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Abstract 

The three-body problem is reexamined in the framework of general relativity. The Newtonian 
three-body problem admits Euler's collinear solution, where three bodies move around the common 
center of mass with the same orbital period and always line up. The solution is unstable. Hence it is 
unlikely that such a simple configuration would exist owing to general relativistic forces dependent 
not only on the masses but also on the velocity of each body. However, we show that the collinear 
solution remains true with a correction to the spatial separation between masses. Relativistic 
corrections to the Sun- Jupiter Lagrange points L\, L2 and L3 are also evaluated. 

PACS numbers: 04.25.Nx, 45.50.Pk, 95.10.Ce, 95.30.Sf 
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Introduction. — The three-body problem in the Newton gravity belongs among classical 
problems in astronomy and physics (e.g, . In 1765, Euler found a collinear solution 

for the restricted three-body problem, where one of three bodies is a test mass. Soon later, 
his solution was extended for a general three-body problem by Lagrange, who also found 
an equilateral triangle solution in 1772. Now, the solutions for the restricted three-body 
problem are called Lagrange points Li,L 2 ,L 3 , L 4 and L 5 , which are described in textbooks 
of classical mechanics [2]. SOHO and WMAP launched by NASA are in operation at the 
Sun-Earth L\ and L 2 , respectively. LISA pathfinder is planned to go to Li . Lagrange points 
have recently attracted renewed interests for relativistic astrophysics js), Q , where they have 
discussed the gravitational radiation reaction on L 4 and L 5 by numerical methods. As a 
pioneering work, Nordtvedt pointed out that the location of the triangular points is very 
sensitive to the ratio of the gravitational mass to the inertial one [5|. Along this course, 
it is interesting as a gravity experiment to discuss the three-body coupling terms at the 
post-Newtonian order, because some of the terms are proportional to a product of three 
masses as M\ x M 2 x M3. Such a term appears only for relativistic three (or more) body 
systems: For a relativistic binary with two masses M\ and M 2 , M\M 2 and M\M 2 exist 
but such three mass products do not. For a Newtonian three-body system, we have only 
the terms proportional to M\M 2 , M 2 M% and M^M\. The relativistic periastron advance of 
the Mercury is detected only after much larger shifts due to Newtonian perturbations by 
other planets such as the Venus and Jupiter are taken into account in the astrometric data 
analysis. In this sense, effects by the three body coupling are worthy to investigate. 

After efforts to find a general solution, Poincare proved that it is impossible to describe 
all the solutions to the three-body problem even for the 1/r potential. Namely, we cannot 
analytically obtain all the solutions. Nevertheless, the number of new solutions is increasing 
jfj]. Therefore, the three-body problem remains unsettled even for Newton gravity. 

The theory of general relativity is currently the most successful gravitational theory 
describing the nature of space and time, and well confirmed by observations. Especially, it 
has passed "classical" tests, such as the deflection of light, the perihelion shift of Mercury 
and the Shapiro time delay, and also a systematic test using the remarkable binary pulsar 
"PSR 1913+16" [7j. It is worthwhile to examine the three-body (or more generally, N- 
body) problem in general relativity. However, it is difficult to work out in general relativity 
compared with Newton gravity, because the Einstein equation is much more complicated [8] 
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(even for a two-body system [9 12] ) . So far, most of post-Newtonian works have focused 
on either compact binaries for an application to gravitational waves astronomy or N-body 
equation of motion (and coordinate systems) in the weak field such as the solar system (e.g. 



13J). In addition, future space astrometric missions such as SIM and GAIA 14H16I] require 



a general relativistic modeling of the solar system within the accuracy of a micro arc-second 



171 ]. Furthermore, a binary plus a third body have been discussed also for perturbations of 
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-|21]. 



gravitational waves induced by the third body 

The Newtonian three-body problem admits Euler's collinear solution, where three bodies 
move around the common center of mass with the same orbital period and always line up. 
The solution is unstable against small displacements. Hence it is unlikely that such a simple 
configuration would exist owing to general relativistic forces dependent not only on the 
masses but also on the velocity of each body. The line could bend at a certain location 
of one mass, which means a V-shape configuration. The above Newtonian instability does 
not necessarily come from small perturbations of acceleration. Therefore, it is interesting 
to ask whether the general relativistic gravity in the rather complicated form admits a 
collinear solution or lead to such a V-shape solution. We shall also evaluate for the first 
time relativistic corrections to L±, L 2 and L 3 for the Sun- Jupiter system. 

In recent, a choreographic solution has been studied in the framework of general relativity 



221 ] . Here, a solution is called choreographic in the celestial mechanics, if every massive 



particles move periodically in a single closed orbit. As a choreographic solution, the figure- 
eight one was found first by Moore and rediscovered with its existence proof by Chenciner 
and Montgomery 23M27I] . The solution was shown to remain true at the first post-Newtonian 



221 ] and also the second post-Newtonian orders [28]. Such an unexpected feature may be 



found in the collinear solution. 

This paper is organized as follows. First, we briefly summarize a usual treatment of 
Euler's collinear solution. Next, we extend the formulation to the post-Newtonian case by 
treating the Einstein- Infeld- Hoffman equation of motion. We take the units of G — c — 1. 
Newtonian Euler's collinear solution. — Let us briefly summarize the derivation of the Euler's 
collinear solution for the circular three-body problem in Newton gravity. We consider Euler's 
solution, for which each mass moves around their common center of mass denoted as Xq with 
a constant angular velocity u. Hence, it is convenient to use the corotating frame with the 
same angular velocity u. We choose an orbital plane normal to the total angular momentum 
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as the x — y plane in such a corotating frame. We locate all the three bodies along a single 
line, along which we take the x-coordinate. The location of each mass Mj (I = 1, 2, 3) is 
written as X I = (x/,0). Without loss of generality, we assume x 3 < x 2 < x%. Let -Re- 
define the relative position of each mass with respective to the center of mass Jg = (xq, 0), 
namely Rj = xi — xq [Ri ^ \X t \ unless xq = 0). We choose x = between M x and M 3 . 
We thus have R 3 < R 2 < R\, R3 < and Rt > 0. 

It is convenient to define a ratio as R23/R12 = z, which is an important variable in the 
following formulation. Then we have R13 = (1 + z)R\2- The equation of motion becomes 



n > M 2 M 3 
jx 12 rt 13 



9 Mi M 3 

rL 12 -^23 



where we define 



o 2 M, M 2 

•"-13 -"-23 

R U = Xj - Xj, (4) 
Ru = \Ru\. (5) 

Figure [1] shows a classical configuration at t — 0. At i = Tjv/2, this configuration is rotated 
by 7r radian, where Tjv denotes the Newtonian orbital period. 

First, we subtract Eq. 02]) from Eq. ([TJ and Eq. (J3J from Eq. fl2]) and use i?i 2 = \Xx— X 2 \ 
and i?23 = \X 2 — X% \. Such a subtraction procedure will be useful also at the post-Newtonian 
order, because we can avoid directly using the post-Newtonian center of mass 



29[ . Next > 

we compute a ratio between them to delete u 2 . Hence we obtain a fifth-order equation as 

(M 1 +M 2 )z 5 + (3M 1 +2M 2 )z 4 + (3M 1 +M 2 )z 3 -(M 2 +3M 3 )z 2 -(2M 2 +3M 3 )z-(M 2 +M 3 ) = 0. 

(6) 

Now we have a condition as z > 0. Descartes' rule of signs (e.g., [3^) states that the number 
of positive roots either equals to that of sign changes in coefficients of a polynomial or less 
than it by a multiple of two. According to this rule, Eq. ([6]) has the only positive root 
z > 0, though such a fifth-order equation cannot be solved in algebraic manners as shown 
by Galois (e.g., [30]). After obtaining z, one can substitute it into a difference, for instance 
between Eqs. (pQ) and (j3J). Hence we get co. 
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FIG. 1: (color online). Schematic figure for a classical configuration of three masses denoted by 
Mi (red), M2 (green) and M3 (blue), which represent a Newtonian collinear solution. The filled 
disks denote each mass at t = 0. Definitions of a, Ri and R[j are also indicated. 



post-Newtonian collinear solution. — In the previous part, the motion of massive bodies 
follows the Newtonian equation of motion. In order to include the dominant part of general 
relativistic effects, we take account of the terms at the first post-Newtonian order. Namely, 
the massive bodies obey the Einstein- Infeld- Hoffman (EIH) equation of motion as Q, 
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where vj denotes the velocity of each mass in an inertial frame and we define 



Ru 



(8) 



Here, the middle term with vector (v A — v K ) has a zero coefficient for the circular collinear 
case, while the remaining accelerations are radial. 

We obtain a lengthy form of the equation of motion for each body. By subtracting the 
post-Newtonian equation of motion for M 3 from that for M\ for instance, we obtain the 
equation as 

R 13 u 2 = F N + F M + F v u 2 , (9) 

where we denote a = Ri 3 and the Newtonian term F N and the post-Newtonian parts F M 
(dependent on the masses only) and Fy (velocity-dependent part divided by u 2 ) are defined 

as 

M 
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respectively. Here, we define the mass ratio as vj = MjjM for M = ^ 7 Mj and frequently 
use 1/2 — 1 — v\ — 1/3. It should be noted that in this truncated calculation we ignore the 
second post-Newtonian (or higher order) contributions so that we can replace, for instance, 
V\ by Riuj (with using the Newtonian R ± ) in post-Newtonian velocity- dependent terms such 
as v\. 

After straightforward but lengthy calculations, which are similar to the above Newtonian 
case, we obtain a seventh-order equation as 

7 

F(z) = J2 A * zk = > ( 13 ) 

fc=0 
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where we define 
A 7 = 

A« = 
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(14) 



(15) 



(16) 



(17) 



(18) 



(19) 



(20) 



(21) 



This seventh-order equation is symmetric for exchanges between v x and u 3 , only if one makes 
a change as z — > 1/z. This symmetry may validate the complicated form of each coefficient. 

Once a positive root for Eq. ([13"]) is found, the root z can be substituted into Eq. (jUJ) in 
order to obtain the angular velocity u. 
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The angular velocity including the post-Newtonian effects is obtained from Eq. (Q as 

" = "» ( 1 + l| + 2^)' < 22 > 
where = (-F/v/i? 13 ) 1//2 denotes the angular velocity of the Newtonian collinear orbit. 

Figure [2] shows a numerical example for Mi : M2 : M3 = 1:2:3, -R12 = 1 and 
a/M = 100, where the post-Newtonian correction is of the order of one percent. In this 
figure, we employ the inertial frame (x, y) but not the corotating frame (x,y). We assume 
23 < X2 < Xi throughout this paper. This figure suggests that as an alternative initial 
condition we can assume x\ < x<i < £3, which is realized at t = T/2 (T=orbital period) in 
this figure. It is natural that this is a consequence of the parity symmetry in our formulation. 
Numerical calculations for this figure show that the relativistic correction in Eq. ()22l) is 
negative, that is u < ujn- It should be noted also that the location of each mass at t = T/2 
is advanced compared with that at t = T/v/2 (a half of the Newtonian orbital period). This 
may correspond to the periastron advance (in circular orbits). 

We produce this figure in two ways. One is that we use our formulation to determine u 
and consequently T. Also un and T/v are obtained at the Newtonian level. Next we rotate 
the configuration by angles ir x (T/v/T) and 7r, respectively. The other is that we directly see 
the evolution of the post-Newtonian system. That is, we solve numerically the EIH equation 
of motion until t = T/v/2 and t = T/2, respectively. The both methods provide the same 
plot. This agreement may also validate our formulation. 

Finally, we focus on the restricted three-body problem so that we can put z = z^(l + e) 
for the Newtonian root zn- Substitution of this into Eq. ( fl3|) gives the post-Newtonian 
correction as 

_ Zfc A PNkZ k N . s 

~ ~V kA, T ,7 k ' [ ' 

2—tk ^^NkZjy 

where A Nk and A PNk denote the Newtonian and post-Newtonian parts of A k , respectively. 
For a binary system of comparable mass stars, the correction e is 0(M/a). This implies 
that a corrected length is of the order of the Schwarzschild radius. 

For the Sun- Jupiter system, general relativistic corrections to L 1; L 2 and L 3 become +30, 
—38, +1 [m], respectively, where the positive sign is chosen along the direction from the Sun 
to the Jupiter. Such corrections suggest a potential role of the general relativistic three (or 
more) body dynamics for high precision astrometry in our solar system and perhaps also 
for gravitational waves astronomy. They are very small but may be marginally within the 
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FIG. 2: (color online). Orbit of each mass representing the post-Newtonian collinear solution in 
an inertial frame. We assume Mi (red) : M2(green) : Ms(blue) = 1:2:3, R12 = 1 and a/M = 100. 
The filled disks denote each mass at t = 0, the triangles at t = Tjv/2, and the circles at t = T/2, 
which can be obtained also by the reflection of the filled disks with respect to the x = line. 

limits of the current technology, since the Lunar Laser Ranging experiment has successfully 
measured the increasing distance of the Moon ~ 3.8cm/yr. 

Conclusion.— We obtained a general relativistic version of Euler's collinear solution for 
the three-body problem at the post-Newtonian order. Studying global properties of the 
seventh-order equation that we have derived is left as future work. 

It is interesting also to include higher post-Newtonian corrections, especially 2.5PN ef- 
fects in order to elucidate the secular evolution of the orbit due to the gravitational radiation 
reaction at the 2.5PN order. One might see probably a shrinking collinear orbit as a conse- 
quence of a decrease in the total energy and angular momentum, if such a radiation reaction 
effect is included. This is a testable prediction. 

It may be important also to search other solutions, notably a relativistic counterpart of 
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the Lagrange's triangle solution (so-called L4 and L5 in the restricted three-body problem). 
Clearly it seems much more complicated to obtain relativistic corrections to the Lagrange 
orbit. 

We are grateful to A. Chenciner, Y. Kojima, K. Tanikawa for useful conversations. This 
work was supported in part (H.A.) by a Japanese Grant-in-Aid for Scientific Research from 
the Ministry of Education, No. 21540252. 
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